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Abstract - Kinetics of homogeneous nucleation and growth of copper precipitates under 
electron irradiation of Fei_2;Cua; alloys at concentrations x from 0.06 to 0.4 at. % and tem- 
peratures T from 290 to 450° C is studied using the kinetic Monte Carlo (KMC) simulations 
and the statistical theory of nucleation (STN). The conventional assumption about the sim- 
ilarity of mechanisms of precipitation under electron irradiation and during thermal aging 
is adopted. The earlier-developed ab initio model of interactions in Fe-Cu alloys is used for 
both the KMC simulations and the STN calculations. Values of the nucleation barrier Fc 
and the prefactor Jo in the Zeldovich-Volmer relation for the nucleation rate J are calculated 
for a number of concentrations and temperatures. For the dilute alloys with x < 0.2%, the 
STN and the KMC results for the nucleation barrier Fc do virtually coincide with each other, 
which seems to confirm a high reliability of the STN for this problem. The STN calculations 
are also used to estimate the temperature dependencies of concentrations which correspond 
to the homogeneous or the heterogeneous precipitation limit, x^°"^{T) and xJ"^*(T), and both 
dependencies are found to be rather sharp. 

Keywords: precipitation, dilute Fe-Cu alloys, kinetic Monte Carlo simulations, statistical 
theory of nucleation 



1. Introduction 

Studies of copper precipitation in irradiated Fe-Cu-based alloys at temperatures 
T ~ 290 — 300° C and low copper concentrations x between x=0.05-0.06% and 
x=0.2-0.3% (in atomic percents, here and below) attract great attention, first of 
all in connection with the problem of hardening and embrittlement of the nuclear 
reactor pressure vessels (RPVs), see, e. g. [H-[§]. Numerous experiments combined 
with some simulations showed that for the usual conditions of electron irradiation, 
that is, when the radiation- induced clusters of vacancies and self-interstitial atoms 
are not formed and the heterogeneous precipitation on them is basically absent, 
SANS and ATP measurements at x < 0.1% typically do not reveal any copper-rich 
precipitates [H, [2, IB]) even though the alloys remain to be strongly supersaturated. 
The absence of a significant homogeneous copper precipitation in these electron- 
irradiated Fe-Cu-based alloys is commonly explained "by the very low volume 
fraction of probably very small precipitates,... and ... one can infer that the limit of 
observable copper precipitation at 290 °C is located between x=0.09 and 0.08%" 
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or "between x=0.1 and 0.2%" [3]. In particular, cluster dynamics simulations 
of homogeneous precipitation at x=0.088%, T = 290° C described in ^ predict "a 
very low density (5-10~^^ m~^) of very big {R ~ 30 nm) copper precipitates". Under 
neutron or ion irradiations when copper precipitation can occur heterogeneously 
at clusters of vacancies or self-interstitial atoms, such precipitation at T = 290°C 
was observed for x > 0.08 [6, i], but not for x < 0.06% 0, Q. 

Therefore, theoretical studies of copper precipitation in very dilute Fe-Cu-based 
alloys near the concentration or temperature limits of such precipitation, xf°"^{T) 
for the homogeneous nucleation and x^^*{T) for the heterogeneous nucleation, seem 
to be interesting from both fundamental and applied standpoints, even for binary 
Fe-Cu alloys which can be considered as model alloys for the RPV steels P, 0,13, |^. 
Such theoretical studies can be particularly interesting if they can sufficiently reli- 
ably predict the temperature dependencies of these precipitation limits at elevated 
T > 290° C for which experimental estimates are either uncertain or absent. 

A consistent ab initio model for studies of precipitation kinetics in Fe-Cu alloys 
has been recently developed by Soisson and Fu [10], and detailed kinetic Monte 
Carlo (KMC) simulations of precipitation in a number of Fe-Cu alloys at different 
T and x based on this model yield a good agreement with available experimental 
data EH- Therefore, this model and KMC methods developed in fl9] seem 
to be prospective to study the above-discussed problem of copper precipitation in 
dilute Fe-Cu alloys. 

We will also use the statistical theory of nucleation and growth of isolated pre- 
cipitates (to be abbreviated as STN) developed by Dobretsov and Vaks [1^. [l^. 
This theory provides microscopic expressions for both the nucleation barrier (ac- 
tivation barrier for the formation of a critical embryo) and the prefactor Jq in 
the phenomenological Zeldovich and Volmer formula for the nucleation rate J (the 
number of supercritical embryos formed in unit volume per unit time) [3] valid 
for the initial, steady-state nucleation stage of precipitation: 

J = Joexp(-/3Fe) (1) 

where j3 = 1/T is the reciprocal temperature. The STN corresponds to a number 
of refinements of earlier models of nucleation suggested by Cahn-Hilliard [la] and 
Langer [l^. It is based on the generalized Gibbs distribution approach [13] and 
enables one to quantitatively calculate values of Fc and Jq for the microscopic 
alloy model chosen, particularly in the case of high nucleation barriers, jiFc S> 1, 
which just corresponds to the dilute alloys under consideration. In description of 
growth of supercritical embryos we also use the phenomenological equation of such 
growth suggested in the classical theory of nucleation [13] and access validity of 
this equation by comparison with our KMC simulations. 

In all our simulations and calculations we consider only homogeneous precipita- 
tion and accept the conventional assumption (discussed, in particular, by Mathon 
et al. 0]) that for the moderate irradiation intensity (i.e. when the ratio of ballistic 
jump frequency to thermally activated jump frequency is small), the main effect of 
irradiation is to increase the point defect concentrations, and hence the Fe and Cu 
diffusion coefficients. It means that the precipitation kinetics under electron irra- 
diation is the same as that during thermal aging at the same temperature, except 
for an acceleration factor Ai^r defined by the relation: 

A,rr = D^SUDcu (2) 

where -D^^ and -D^^ are the copper diffusivities during thermal aging and under 



7 Philosophical Magazine VSZ-final 



Philosophical Magazine 3 

irradiation, respectively. As evolution of microstructure under precipitation is de- 
termined by the diffusion of copper, Eq. ([2]) implies an analogous relation between 
the evolution times tirr and tth under irradiation and during thermal aging: 

tirr — tifi/Aij-j-. (3) 

Then the precipitation kinetics under electron irradiation can be described by the 



KMC codes developed by Soisson and Fu [lO(] for thermal aging with replacing the 



thermal aging time tth by the "scaled" time tirr according to Eq. ([3]) . This approach 
is accepted throughout this work. In more detail, scaling relation ([2]) and estimates 
of the acceleration factor Airr are discussed in Sec. 4. At the same time, the main 
results and conclusions of this work do not depend on A^rr values as this scaling 
factor is canceled in many important kinetic characteristics. 

In Sec. 2 we present some results of calculations of structure of critical embryos 
and parameters Fc and Jq in ([T]) for dilute Fe-Cu alloys obtained using the STN 



mm 

In Sec. 3 we describe our KMC simulations of nucleation and growth of 



copper precipitates in dilute Fe-Cu alloys. In Sec. 4 we present an estimate of the 



acceleration factor A^rr based on the rate theory models [221 ] and find this estimate 
to reasonably agree with those obtained using the available experimental data [2] 
and our KMC simulations. In Sec. 5 we use the phenomenological equation for 
growth of a supercritical embryo and our KMC simulations to estimate some 
important kinetic parameter which enters into the STN expression for the prefactor 
Jo in Eq. ([1]). In Sec. 6 we show that the nucleation barrier Fc values calculated us- 
ing the STN agree well with those estimated in our KMC simulations, particularly 
at large /3Fc ^ 1 which correspond to the concentrations and temperatures near 
the above-mentioned precipitation limits. In Sec. 7 we use the STN calculations 
combined with some plausible physical assumptions to estimate the temperature 
dependence of both the homogeneous and the heterogeneous precipitation limit, 
x^°"^{T) and a;^^*(T), for practically interesting temperatures T between 290 and 
450°C. Our main conclusions are summarized in Sec. 8. 



2. Properties of critical embryos in dilute Fe-Cu alloys calculated using the 
statistical theory of nucleation 

In this section we use the STN described in Refs. [13, [13] (to be referred to as I 
and II) to calculate some thermodynamic, structural and kinetic characteristics of 
critical embryos. For these calculations we use the ab initio model of Fe-Cu alloys 
developed by Soisson and Fu and described in detail in Ref. [10]. Here we only note 
that this model uses the following values of the binding energy between two copper 
atoms and between a copper atom and a vacancy, -EcuCu ^^'^ ^vCu^ n-ih 
neighbors (in eV): 

E^^^Cu = 0.121 - 0.182r, E^J^cu = 0-021 - 0.09ir, 

Eth, = 0.126, Etl, = 0.139. (4) 

The high values of Eq^q^ in Eqs. ([4j) correspond to strong thermodynamic driving 
forces for precipitation, while a strong attraction between vacancy and copper 
atoms leads to the strong trapping of vacancies by copper precipitates discussed in 
detail in [lo| . 

First we consider the structural and thermodynamic characteristics of embryos 
which are described in the STN in terms of the occupation number rii for each 



:7 Philosophical Magazine VSZ-final 



4 Taylor & Francis and I. T. Consultant 

lattice site i. The operator rii is unity when a copper atom is at site i and zero 
otherwise, while probabilities of various distributions {n^} are described by the 
distribution function P{nj} given by Eq. (1-3): 

P{ni} = exp[/3in + fiY,ni-H)]. (5) 

i 

Here H is the configurational Hamiltonian supposed to be pairwise, ^ is the chem- 
ical potential, and the grand canonical potential O is determined by the normal- 
ization condition: 

H = VijUiUj, = — TlnTr exp[/3(^^ A^nj — H)] (6) 

i>j i 

where Vij are the configurational interactions, while symbol "Tr" means summa- 
tion over all configurations {n,}. The free energy F of an alloy with embryo is 
determined by relations 

F = + 11 Q , /X = dF/dci = constant (7) 

i 

where Cj = (n,) = Tr(njP) is the mean occupation number, or the local concen- 
tration. The free energy F = F{ci} in ([7|) was calculated using several statistical 
approximations described in detail in I and II: mean-field, "mean-field with fluc- 
tuations", pair-cluster, and "pair-cluster- with fluctuations". The latter was found 
to be more consistent and accurate than others, thus below we use only the "pair- 
cluster-with fluctuations" approximation. 

In tables 1 and 2 and Fig. 1 we present some characteristics of critical embryos. 
To illustrate both concentration and temperature dependencies of these character- 
istics, we consider two series of Fei-xCu^: alloys: those with the same temperature 
T = 290° C but different concentrations x varying from 0.06 to 0.4% (Table 1), and 
those with the same concentration x = 0.2% but different temperatures T vary- 
ing from 290 to 390° C (Table 2). In these tables, Nc is the total number of copper 
atoms within the embryo, and Rc is its effective radius defined as that of the sphere 
having the same volume as Nc copper atoms in the BCC lattice of a-iron with the 
lattice constant ao=0.288 nm: 

= ao(3A^c/8vr)^/=^ = 0.142 N}'^ nm. (8) 

Note that this critical radius (as well as the precipitate radius R in Eq. ()27p be- 
low) characterizes the total number of atoms but not the "geometrical" size of 
precipitate, in particular, not its mean squared radius {rf)^/"^ defined by Eq. (II- 
10). Therefore, for small and "loose" precipitates, such as those shown by lines B 
and C in Fig. 1, these radii can notably differ from "geometrical" ones. However, 
quantities Rc and i? in ([8]) and ([27l) are convenient to describe the precipitation 
kinetics. 

Quantities AJIq and Af^i in tables 1 and 2 are the zero-order and the first-order 
terms in fluctuative contributions to the total nucleation barrier Fc = Af^o -|- Al^i 
(defined as the difference between grand canonical potentials of two alloy states, 
that with the embryo and that without the embryo). These two terms correspond 
to the iterative treatment of some exact relation of thermodynamics of nonuniform 
systems given by Eq. (1-9) which relates the free energy F = F{ci} to the correlator 
Kij of fluctuations of site occupations ni: Kij = {{rii — Ci){nj — cj)). As discussed 
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in II and below, the accuracy of such iterative treatment of fluctuations for dilute 
alloys considered is usually rather high. Fig. 1 illustrates variations of structure of 
critical embryos with concentration or temperature. 

The results presented in tables 1, 2 and Fig.l clearly illustrate very sharp varia- 
tions of both structure and thermodynamics of critical embryos with concentration 
X and temperature T. In particular, at T = 290°C, the decrease of concentration 
from X = 0.4 to a; = 0.06% leads to the increase of the embryo size Nc and the 
reduced nucleation barrier /3Fc by about three times: from Nc — 8 to Nc ~ 25, 
and from f3Fc — 12 to (3Fc — 35. According to the Zeldovich-Volmer relation ([T|), 
the latter implies decreasing the nucleation rate by about ten orders of magnitude 
(variation with x of the prefactor Jq in ([T|) will be shown to be negligible). Table 
2 illustrates a similar sharp increase of Nc and /3Fc under elevating temperature 
between 290 and 390°C. Fig. 1 shows that this sharp rise of sizes and nucleation 
barriers is accompanied by notable changes in the structure of the embryo: its 
boundary at high PFc > 30 is much less diffuse than that at moderate f3Fc < 20. 

Discussing the concentration profiles shown in Fig. 1 we note that these profiles 



(as well as those shown in Fig. 1 in [ll|] and Fig. 2 in [18]) usually correspond to 
rather diffuse interfaces. At the same time, KMC simulations show that only very 
small clusters can have a diffuse interface, while for R > 0.3 nm (that is, A'^ > 10 
copper atoms) one observes almost "pure" copper clusters with a sharp interface, 
see, e.g.. Fig. 8 in [10]. However, one should take into account that the profiles 
shown in Fig. 1 (and other similar figures) correspond to statistical averaging over 
all orientations of a cluster which is typically rather anisotropic, and this averaging 
leads to diffuse interfaces; it is illustrated, in particular, by Fig. 8 in [l3]. Let us also 
note that in phenomenological treatments of critical embryos based on the Cahn- 
Hilliard continuum approach jl^, 20], the embryo interfaces are usually notably 
more diffuse than those obtained in the microscopic STN used. In particular, rather 
sharp interfaces shown by curves A and D in Fig. 1 can hardly be obtained in the 
phenomenological treatments. 

Let us now discuss the prefactor Jq in Eq. ([T]). In the STN, it is given by Eq. 
(n-3): 

Jo = (/3l7ol/2vr)i/'A/'L»R(u)Da,. (9) 

Here the first three factors (discussed in detail in I) have the "thermodynamic" 
origin. They describe dependencies of the free energy and the distribution function 
of the embryo on the "critical" variable a characterizing its size, and on three 
variables u characterizing tits position R. The last factor Daa is the generalized 
diffusivity in the Q-space that corresponds to an increase of a, i. e. to growth of the 
embryo. This factor is defined by Eq. (1-55) which expresses Daa as a certain linear 
combination of generalized mobilities Mij which describe the temporal evolution 
of mean occupations Cj via Eq. (1-50): 

dc^/dt = [Mij - 6i, J2 Mik]f3dF/dcj . (10) 



The Mij values can be calculated using some microscopic models, for example, 
those employed in II. In particular, for the ideal solution (or in the dilute alloy 
limit which corresponds to small q), the mobility Mij = Mf- is given by Eq. (II-8) 
with itjfe = 0: 

= %,nn7LN(l - Ci)c,{l - C,)fl\ (11) 
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Here 5ij^nn is the Kroneker symbol equal to unity when sites i and j are the nearest 
neighbors and zero otherwise, while quantity 7^^ has the meaning of the mean 
rate of exchanges between neighboring copper and iron atoms. It is related to the 
copper diffusivity Dq in a pure iron by the following relation (see, e. g., Eqs. (72) 
and (74) in ^23;]): 

7^ = Do/al (12) 

where ao is the BCC iron lattice constant used in Eq. ([8]). 

At the same time, mobilities Mij can be considered as phenomenological param- 
eters in the general Onsager-type equations (jlOp . and for weakly nonuniform alloys 
under consideration they can be supposed to be nonzero only for neighboring sites 
i and j, just as in Eq. (|lip . For our problem, the phenomenological model of Mij 
can be considered as adequate if it properly describes the diffusion of solute (cop- 
per) atoms in processes of nucleation and growth of embryos. In Sec. 5 we show 
that growth of a supercritical embryo in dilute Fe-Cu alloys can be described well 
using some generalization of the phe nomenological diffusion equation suggested in 
the classical theory of nucleation [ij] given by Eq. (j26p below. This generalized dif- 
fusion equation differs from those for ideal or dilute solutions mainly by the value 
of the effective diffusivity Df,ff which takes into account the acceleration of pre- 
cipitation kinetics due to the high mobility of small copper clusters noted in Refs. 



I21I ] and discussed below. The i^e// values for dilute Fe-Cu alloys considered 
are estimated in Sec. 5 using kinetic Monte Carlo simulations. 

Therefore, we suggest that the adequate phenomenological description of pre- 
cipitation in Fe-Cu alloys under consideration can be obtained if we assume for 
mobilities Mij in (llOp the same simplest form as that in Eq. (jlip . 

Mij = 6ij^nnlnn[Ci{l - Ci)Cj{l - Cj)]^/^ , (13) 

but assume the effective exchange rate 7^™ in (llSp to be proportional not to the 
diffusivity Dq in a pure iron as in Eq. (I12p . but to the phenomenological effective 
diffusivity -De// mentioned above: 

Inn = Deff/al. (14) 

Validity of this assumption will be checked by comparison with the KMC simulation 
results presented in Sees. 5 and 6. 

In Table 3 we present values of various factors in the expression Q for the pref- 
actor Jo in Eq. ([T]) for some alloys Fe-Cu. Quantity 70 in this table is the derivative 
of the free energy F with respect to the effective embryo size a: 70 = d'^F/do?. 
As the critical embryo corresponds to the saddle-point of F in the Q-space with 
respect to this size (I, II), 70 values in table 3 are negative. M is the normalizing 
constant in the embryo size distribution function; Daa is the generalized diffusivity 
discussed above, and /^^(u) is some geometrical factor which for a large embryo 
is proportional to its surface. For quantities Daa-, -Dr(u), and the total prefac- 
tor Jo we present their "reduced", dimensionless values: Daa /inn, i'a-CR(u), and 
■-^0 = VaJo/lnn- Note that the embryo size Nc decreases with x and increases with 
T, as tables 1 and 2 show. 

Table 3 shows that both concentration and temperature dependencies of quan- 
tities TV and Daa/lnn are rather weak. On the contrary, the geometrical factor 
Dr,(u) notably increases with the embryo size while /3|7o| somewhat decreases 
with Nc- The total reduced prefactor Jq = VaJ^/inn notably increases with Nc, by 
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about 20 or 5 times for the total x or T intervals shown in table 3. However, in 
the total nucleation rate J given by Eq. ([T]), these variations of Jq are negligible as 
compared to the above-mentioned huge changes of the activation factor exp {—/3Fc) 
for these intervals of x or T. 

Now we note that a very sharp decrease of this activation factor with x be- 
tween X ~ 0.2 and x ~ 0.08% seen in table 1 evidently correlates with the 
above-mentioned experimental estimates of the homogeneous precipitation limit 
at T = 290° C given in Refs. UMM- 

It seems natural to suggest that the position 
^hom q£ ^j^jg limit is mainly determined by the activation factor exp(— /3Fc), or 
by the reduced nucleation barrier l3Fc which at x = x^°^ reaches a certain high 
value {PFc)i. For example, if we accept for the experimental x^°"^=xi at 290°C 
the estimate from Ref. [2], x^^^ ~ 0.08-0.09%, we have: {l3Fc)i ~ 28, while taking 
this estimate from Ref. p], rr^^^ ~ 0.15%, we obtain: {j3Fc)i — 23. In more detail, 
estimates of the precipitation limits xi are discussed below in Sec. 6. 



3. Kinetic Monte Carlo simulations of nucleation and growth of precipitates 
in dilute Fe-Cu alloys 

As mentioned in Sec. 1, our KMC simulations use the "scaling" assumptions ([2]) 
and ([3]) which enable us to relate the homogeneous precipitation in irradiated alloys 
to that during thermal aging for which the KMC codes developed by Soisson and 
Pu can be employed. In these simulations we consider the following Fei-^Cu^; 
alloys: (A) those at T = 290° C and x equal to 0.088, 0.116, 0.144, 0.177, 0.2, 0.3 
and 0.4%, and (B) those at x = 0.2% and temperatures T equal to 290, 305, 320, 
335, 350, 365, 380 and 390° C. We use simulation volume Vg = Lp' with periodic 
boundary conditions and, to check the statistical reliability of results, we usually 
employ several different values of the simulation size L: 

Li = 128oo, L2 = 160ao, L3 = 192ao. (15) 

Before to describe the results, we mention a characteristic feature of such simu- 
lations for some low concentrations x or high temperatures T: at such x and T, we 
observe only one precipitate within simulation volume, that is, the total number 
Np of precipitates within simulation box is unity (while at further lowering x or 
elevating T we observe no precipitates for the simulation time). In our simulations 
it was the case for the series (A) alloys with T = 290°C at x < 0.144%, and for the 
series (B) alloys with x = 0.2%, at T > 335°C. In all such cases, the relation Np = 1 
preserves under increasing simulation volume (replacing Li by L2 or L3), though 
the incubation time tine (that which precedes to the formation of the embryo [2j]) 
somewhat decreases when Vg increases. 

The presence of changes of physical characteristics (such as the precipitate den- 
sity dp = Np/Vs or the incubation time tj„c) under variations of simulation volume 
evidently indicates on the statistical unreliability of simulation results for these 
characteristics. At the same time, the growth of a precipitate after its formation 
seems to be described by such simulations quite properly. Therefore, even though 
the simulations with Np = 1 can not be used to study the evolution of precipitate 
density dp{t), they will be widely used for studies of precipitate growth described 
in Sec. 5. 

When the number of precipitates within simulation box significantly exceeds 
unity, the simulations become statistically reliable and can be used to study the 
precipitate density dp{t). In figures 2 and 3 we show the dependencies dp{t) for all 
our simulations for which the final number of precipitates within simulation box is 
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not too small: Np > 3. The simulation time on abscissa axis is given in the "thermal 
aging" values tth which in our model are related to the observed time tirr by Eq. 
([3]). The acceleration factor Ai^r in this relation and comparison with the available 
data [2] about precipitation in an Fe-0.3%Cu alloy under electron irradiation are 
discussed below in Sec. 4. In the lower part of frames 2(e) and 2(f) we also show 
temporal dependencies of the mean precipitate radius Rm for two simulations with 
X = 0.3%. For other our simulations with Np > 10 these dependencies are similar 
and describe usually a smooth increase of Rm by about 20-30%. 

The dependencies dp{tth) presented in Figs. 2 and 3 are approximately linear. 
Therefore, if we define the nucleation rate J = Jkmc for these simulations as the 
ratio 



where r = (t^/j — t^^'^) is the evolution time counted off the incubation time t^^^, this 
nucleation rate for each simulation remains approximately constant, particularly 
for the simulations with not small Np where fluctuations are not too strong. Hence 
these our simulations seem to describe mainly the first stage of precipitation, the 
steady-state nucleation, for which the constant value of the nucleation rate J is 
characteristic. 

The dashed line in each frame of Figs. 2 and 3 shows our estimate of the Jkmc 
value in Eq. (jl6p assuming it to be constant. For this estimate we tried to use the 
maximum broad interval of time tt^ for which the temporal dependence dp{tth) is 
close to linear. In spite of the evident scatter of KMC results presented in Figs 2 and 
3, particularly for simulations with not large Np, the linear temporal dependence 
dp{Tth) can usually be followed rather clearly . Therefore, our estimates of Jkmc 
values in Eq. (jl6p seem to be sufficiently definite. 

At the same time, the Zeldovich-Volmer relation ([1]) is valid just for the steady- 
state nucleation stage. Therefore, we can use our estimates of Jkmc to estimate 
the parameters Fc and Jq in ([T|) putting: 



As discussed in Sec. 2, variations of the activation factor exp(— /3Fc) with con- 
centration or temperature affect the total nucleation rate J much stronger than 
those of the prefactor Jq. Therefore, for the given JkmC) the activation factor 
exp(— /?Fc) can be determined from Eq. ()17p rather accurately even if the prefactor 
Jo is estimated not too precisely, in particular, when Jq is estimated from the STN 
calculations described in Sec. 2 which include model assumptions ()13p and (I14p . 

Substituting for the prefactor Jq its expression via quantities Jq , Va and "^nn given 
in the last line of table 3, and using Eq. (fTljl which relates ^nn to the effective 
diffusivity -De// mentioned in Sec. 2, we can re- write relation (fT7|) as the equation 
for the reduced nucleation barrier pPc.: 



where index "KMC" in the left-hand side indicates that this expression for is 
based on the KMC simulations. 

Values of the reduced prefactor Jq can be taken from our STN calculations 
illustrated by the last line of table 3. Hence to find the reduced nucleation barrier 
from Eq. (jlSp , we need the effective diffusivity Dg/ / • Estimates of D^f / based 
on our KMC simulations of growth of precipitates are described below in Sec. 5. 



Jkmc = dp{T)/T 



(16) 



Joexp(-/3Fc) = Jkmc- 



(17) 




(18) 
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4. Estimate of acceleration precipitation factor Ai 
of dilute iron-copper alloys 

The concentrations of point defects (vacancies and interstitials) under permanent 
irradiation niay frequently exceed their equilibrium values by several orders of 
magnitudes [2^. In such conditions diffusive phase transformations, such as pre- 
cipitation or ordering, are strongly accelerated. For the case of electron irradiation 
of Fe-Cu alloys this has been observed, in particular, by Mathon et al. [2] and by 
Radiguet et al. @|. If other irradiation effects (such as ballistic mixing or radiation 
induced segregation) can be neglected, it is natural to expect that the precipitation 
kinetics under irradiation is the same as that during thermal aging at the same 
temperature, except for the strong acceleration of this kinetics described by the 
fcictor 74.27-7' in Eqs. ([2]) and ([3]). In this section we estimate this acceleration fac- 
tor for dilute Fe-Cu alloys under typical conditions of electron irradiation [5] and 
compare this estimate with both the experimental and KMC simulation results. 

For dilute alloys under consideration, the copper diffusivity under thermal aging 
in Eq. ([2]) is given by the conventional expression [25| : 



D'(^^ = a,cim,. (19) 

Here is the equilibrium vacancy concentration which is expressed via the en- 
thalpy and entropy of vacancy formation, hI°^ and Sl°^ ^ as: Cv'^ = exp{Si"^ — 
pni"^), and Dy is the vacancy diffusion coefficient in a pure iron which is related 
to the vacancy migration enthalpy and entropy, H^^^ and S^^^ , as follows: 

Dy = aluo exp(5r^ - (20) 

where oq is the BCC lattice constant and uq is the attempt frequency supposed 
to have the order of the Debye frequency. The coefficient ay in (fT9|) describes 
vacancy-solute correlations, and it is commonly written as [2^: 

«. = /2exp(/3Gg^„) (21) 

where /2 is the impurity correlation factor, and G^JJ_^ is the copper-vacancy bind- 
ing energy. These parameters may be obtained from experimental measurements, 
or from ab initio calculations [13] ■ Similarly, under irradiation one may write: 

D'S: = ayct'Dy + aicrDi, (22) 

where d^^ and c-'"'' are the vacancy and the interstitial concentrations under irra- 
diation. 

It seems natural to assume that the copper-vacancy correlation effects and bind- 
ing are not modified by irradiation. Detailed information on the interstitial diffusion 
(needed for calculations of Oi and Di) is usually more difficult to obtain. However, 
when the point defect concentrations reach their steady-state value, one can show 
that dy'^Dy = cf^Di [22,]. If ay and are of the same order of magnitude, one fi- 
nally gets ~ 2aydy'^ Dy, which corresponds to the following acceleration factor 
in di) and Q: 

A,rr = 2{cl''/cl'^). (23) 



To compare our KMC simulations during thermal aging and the experimental 
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kinetics of precipitation under irradiation, the key point is therefore to get a reliable 
estimation of the vacancy poin t defect concentrations under irradiation. According 
to the rate theory models [22], if point defects created by irradiation disappear by 
mutual recombination or by annihilation at dislocations (which are assumed to be 
the dominant point defect sinks), the evolution of their concentrations in a pure 
metal is described by the following equations: 

= K - RciCy - PdDv{cv - cl'^) 

dci 

= K - RciCy - pdDi{ci - cl'^) ~ - RciCv - pdDia . (24) 

Here K is the rate of formation of Frenkel pairs under electron irradiation (in 
dpa/s), pd is the dislocation density, and R is the vacancy-interstitial recombination 
rate: 

R = ATTryi{Dy + Di)/va (25) 

where r„j is the recombination radius, and Va = aQ/2 is the atomic volume. 

For pure iron (for which lattice constant is ao=0.288 nm), most of the point 
defect properties have been estimated by ab initio calculations. In particular, for 
the vacancy formation parameters we have: hI°^ = 2.18 eV [lo| and 5^^°'' = 4.08 
261], which gives: c%'^ ~ 1.8 x 10~^^ at 290°C. The migration enthalpies are: ff™^ = 
0.68 eV [l^ and iJf*^ = 0.34 eV [271]. The Debye frequency is uq = lO^^g-i^ ^nd 
with these parameters a vacancy migration entropy 5™^ ~ 2.2 is required to get 
the experimental pre-exponential factor for self-diffusion of iron (see Refs. 0, m 
for details). Lacking both experimental and ah initio estimations, we will assume 
the same value for the migration entropy of the interstitials. 

The irradiation conditions of ref. are: K = 2 y. 10~^ dpa/s at 290°C, with 
dislocation densities between pd = 10*^ cm^^ and pd = 10^^ cm^^. The evolution 
of point defect concentrations, obtained by numerical integration of Eq. (j24p with 
these two dislocations densities, is given on Fig. 4. In both cases, the steady-state 
values are reached very rapidly (after less than 10~^ s), with vacancy concentrations 
^trr ^ 2,,^ X 10~^^ and 3.6 x 10~^^. The acceleration factor is therefore Ai^r ~ 4 x 10^ 
for the highest dislocation density and Aij-r — 3.9 x 10'' for the lowest. 

In Ref. [2], a precipitate density of 0.9 x lO^^m"^ is observed under electron irra- 
diation of Fe-0.3%Cu alloy after approximately 8.3 x 10^ s. Our KMC simulations 
of precipitation during thermal aging shown in Fig. 2(f) predict a similar density 
after 1.5 x 10^ s. This corresponds to an acceleration by a factor Airr ~ 1.8 x 10^ 
under irradiation. Considering the strong approximations used in the rate theory 
model, and the uncertainties on the point defect properties, the agreement between 
our estimate of Airr and the combination of experimental and KMC results seems 
to be reasonable. 

A more detailed comparison of our simulations with experiments by Mathon et 
al. i is hindered by a low resolution in measurements of precipitate sizes R by 
SANS method used in [3]: Ai? ~ 0.5 nm. It can explain quantitative disagreements 
between the results reported in [3] and those shown in frames 2(e) and 2(f). The 
maximum value of the precipitate density dp{t) reported in [3], dmax ~ 1-65 x 
10^^1/m'^, is by about 1.5 times lower than the final values dp{t) shown in Figs. 2(e) 
and 2(f), even though the maximum value dmax is probably not reached yet in these 
simulations. Hence, our simulated dp exceed those reported in 0] by at least several 
times. At the same time, values of the mean precipitate radius for the nucleation 
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stage (which corresponds to dp{t) < dmax) reported in j^.j ~ 1.1 — 1.65 nm, 

are by about 2-3 times hi ghe r than the analogous values ~ 0.3—0.5 nm observed 
in both experiments and simulations (lol . [ll| for other Fe-Cu alloys, while 



our Rm in frames 2(e) and 2(f) are close to these usual R^. Therefore, it seems 
probable that because of the above-mentioned low resolution Ai? ~ 0.5 nm, values 
of the precipitate density dp reported in (2] are underestimated, while those of the 
mean precipitate size Rm Q-i^s overestimated. More accurate measurements of dp 
and Rm in dilute Fe-Cu alloys under electron irradiation are evidently needed for 
a quantitative comparison with our simulations. 



5. Estimates of effective diffusivity of copper atoms for growth of a 
precipitate using kinetic Monte Carlo simulations 

In the phenomenological theory of growth of a supercritical embryo developed for 
a spherical embryo of a radius R with a sharp edge, the temporal dependence R{t ) 
(for the rigid lattice alloy model used) is described by the following equation [14l ]: 

dR/dt = DeffiR- Rc)ic- Cb)/R^. (26) 

Here -De// is the effective diffusivity of a solute atom (in our case, of a copper atom) 
in an alloy; Rc is the critical radius; c = x is the copper concentration far from 
the embryo; and Cb is the binodal concentration (solubility limit) for temperature 
under consideration. 

In more realistic descriptions, in particular, in the STN or KMC simulations 
described in sections 2 and 3, both critical and supercritical embryos have not 
a sharp edge but a diffuse surface illustrated by Fig. 1. More important, in the 
phenomenological derivation [13] of Eq. ([26]) . the diffusivity D^ff was supposed to 
be constant independent of the precipitate size and its mobility m{N). At the 
same time, the recent KMC studies 10, 21 1 have shown that the small precipitates 



under consideration which contain ^100 copper atoms are highly mobile, being 
much more mobile than individual copper atoms. It leads to a great acceleration of 
precipitation kinetics which qualitatively corresponds to an increase of the effective 
diffusivity D^jf in Eq. ([26]) with respect to the diffusivity of an individual copper 
atom. 

In spite of all these simplifications, the phenomenological equation (j26p is com- 
monly believed to realistically describe growth of supercritical embryos [l^. Ac- 
cepting this point of view, we first present the explicit solution of Eq. (|26p for R{t). 
Then we use our KMC simulations to access validity of Eq. (|26p and to estimate 
the parameter -De// in this equation for dilute Fe-Cu alloys under consideration. 

First we note that employing the simplest sharp-edge model of an embryo in 
usual derivations of Eq. (f26]l seems to be unessential, and this equation (at large 
Nc ^ 1 considered) appears to be valid for any realistic description of growth of 
an embryo, including that used in KMC simulations. The radius R{t) in Eq. (I26p 
mainly characterizes the total number N of solute atoms within the embryo, and 
for any precipitate this radius can be expressed via N analogously to Eq. ([8]) for 
the critical embryo: 

R = ao(3A^/87r)^/3 = 0.142 A^^/^ nm. (27) 

The right-hand-side of Eq. (126p contains two basic factors (in addition to the dif- 
fusivity -De// and a geometrical factor 1/R'^) which naturally describe the driving 
force for growth of the embryo: the factor (c—Ch) that characterizes supersaturation 
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of a metastable alloy, and the factor (R — Rc) that describes vanishing of this force 
at i? = Rc- Therefore, the phenomenological equation (j26l) with the generalization 
(I27p seems to be rather plausible, and comparison to our KMC simulations given 
below seems to confirm its validity, at least for the R and t intervals studied. 
First we present the explicit solution R{t) of Eq. (f26]) . As discussed in detail in 



IJ], the embryo can be considered as "supercritical" only when its size R exceeds 
some value {Rc + A) where A determines the scale of critical fluctuations of sizes 
near R^.- Therefore, the time tj when the supercritical embryo has been eventually 
formed is defined by the relation: 

= (i?c + A). (28) 

For the sharp-edge model of the embryo used in [14], the fluctuation width A 
is expressed via the interfacial energy cr as: A = (T/Svrcr)^/^. To generalize this 
estimate to the case of real embryos with diffuse interfaces, we can express a via 
the critical radius Rc and the nucleation barrier Fc- It yields: 

A = Rca, a = {IT/Wcf^. (29) 

For large embryos under consideration with (3Fc ^ 1, values of a are small; in 
particular, for the Fe-Cu alloys listed in tables 1 and 2, we have: a ~ 0.15 — 0.2. 

Integrating Eq. (|26p from the initial time t f defined by Eq. (j28p to the arbitrary 
time t when R{t) = R, we obtain: 

{R-Rc- A){R + 3Rc + A)/2 + R^ In [{R - i?,)/A] = Dcff{c - Cb)r (30) 

where r = {t — tf) is the total time of growth of a supercritical embryo. If we 
describe this growth by a reduced variable y = {R — Rc — A)/Rc, Eq. ([30]) takes a 
universal form which contains only the reduced dimensionless time 

y(2 + a + y/2) + ln(l + y/a) = 

C = Dcff{c-ct)T/Rl (31) 
At large R 3> Rc, the dependence R{t) given by Eq. ([30]) or (j3T|) takes the form 

RiT) = [2Dcff{c-c,)T]y^ (32) 
which describes the diffusion-controlled growth of a spherical embryo, see, e. g.. 



Now we compare the description of growth of an embryo by the phenomenological 
equation (I30p or ()3ip with that given by our KMC simulations. First, we note that 
for our model of irradiated alloys based on Eqs. ([2]) and ([3]), the reduced time ^ in 
Eq. (I3ip does not include the irradiation acceleration factor Airr, just as the right- 
hand side of Eq. (jlSp , and this reduced time is explicitly expressed via the thermal 
aging time t = Tth = {tth — i{h) used in our KMC simulations. Therefore, growth 
of a supercritical embryo can be described, from one side, by the function R{Cj 
determined by Eq. ([30]) or ([3T]) . From the other side, such growth can be followed 
in our KMC simulations for any precipitate chosen, which yields the dependence 
R = RKMc{tth) with R defined by Eq. ([27]) . Equating these two quantities. 



= RKMcitth)i 



(33) 
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with ^ given by Eq. (I3ip for r = Tth, we obtain the equahty which contains only 
one unknown parameter -De// at aU times Tth > 0. 

One can expect that the proper choice of this single parameter can provide a good 
accuracy for obeying equality (I33p at all Tth considered only if the phenomenological 
equation (I26p holds true for these Tth- Therefore, the check of validity of relation 
(p3]) enables us, first, to access reliability of the phenomenological equations (f26]l 
and (|30l) for description of growth of supercritical embryos and, second, to estimate 
the effective diffusivity -De// in these equations for dilute Fe-Cu alloys. 

It is convenient to characterize the effective diffusivity Deff{x,T) by its ratio Ac 
to the appropriate dilute alloy diffusivity D{0,T) = Dq{T) writing D^jf as 

Dcff{x,T) = AcXDoiT). (34) 

The parameter Ac = Ac{x,T) in (I34p characterizes the effective acceleration of 
growth of precipitates being mainly due to the above-mentioned high mobility 
of small copper clusters. To differ it from the irradiation acceleration factor Airr 
discussed in Sec. 4, Ac will be called "the acceleration diffusion parameter". 

For the dilute alloy diffusivity Dq(T) of a copper atom in a pure iron at T = 
290 — 450° C we will use the conventional Arrhenius-type expression 

Do{T) = Aexp{-PQ) (35) 

with the A and Q values suggested by Soisson and Fu [l^ basing on combinations 
of their ab initio calculations and empirical estimates: 

A = 0.97 • IQ-^ mVs , Q = 2.67 eV. (36) 

In Figs. 5 and 6 we show the phenomenological functions R{^) defined by Eqs. 
(pOj) and (f3T]l together with the dependencies RKMcitth) observed in our KMC 
simulations. Each curve RKMciUh) in these figures starts from the time of formation 
of the precipitate chosen observed in the KMC simulation, while the curve R{S,) = 
R[^{tth)] starts from the time of formation of the supercritical embryo defined 
by Eq. (|28p . Each frame in these figures corresponds to some of our simulations; 
for other simulations (not shown in Figs. 5 and 6), the results are similar. The 
parameter Ac and the time used to draw the curve R{^) in each frame have been 
estimated from the best fit of this R{(,) to the RKMcitth) curve at this frame. In 
Fig. 7 we show the concentration and temperature dependencies of the acceleration 
diffusion parameter Ac obtained in these KMC estimates. 

Let us discuss the results presented in Figs. 5-7. First, we note that some of 
irregularities in dependencies RKMc{tth) seen in Figs. 5 and 6 can be related not to 
the real statistical fluctuations of sizes but to the methodical errors due to too long 
temporal intervals between sequent savings of KMC data. This seems to be the case, 
in particular, for initial stages of simulations shown in frames 5(a), 6(c) and 6(d), 
and also for some other frames, e. g., 6(a). At the same time, for the later stages of 
growth when the precipitate becomes supercritical, R > (Rc + A), these methodical 
distortions seem to be less significant. Second, Figs. 5 and 6 clearly illustrate the 
strong fluctuations of sizes in the "critical" region R < (Rc + A) mentioned in the 
discussion of Eqs. ([28]) and (f29l) . Such fluctuations are pronounced, in particular, 
in frames 5(b), 5(c), 6(c), 6(e), and 6(f). Third and most important. Figs. 5 and 
6 show that the description of growth of precipitates by the phenomenological 
equations (I26p and (I30p seems to agree well with the KMC simulations, at least up 
to R> 1.5Rc, that is, under increase of the number of copper atoms in the growing 
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embryo by about four times. 

Figs. 5 and 6 also clearly illustrate the effects of cluster mobility and their direct 
coagulation mentioned above. In particular, sudden increases ("jumps") of R{t) 
seen in Fig. 5(d) at tth = 10.6 x 10^ s, as well as in Fig. 5(f) at tth = 1.5 x 10^ and 
tth = 1.75 X 10^ s, occur simultaneously with the disappearance (evidently, due to 
coagulation) of one of supercritical clusters in Figs. 2(b) and 2(f). Similarly, jumps 
of R{t) in Fig. 6(a) at tth = 6.56 x 10^ s and in Fig. 6(b) at tth = 1-46 x 10^ s 
occur simultaneously with the disappearance of one supercritical cluster in Figs. 
3(b) and 3(d), respectively. The smaller jumps in R{t) seen in Figs. 5(d), 5(e), 
6(a) and 6(b) can correspond to coagulation of subcritical clusters which are not 
registered in Figs. 2 and 3. At the same time, for simulations with Np = 1 shown 
in Figs. 5(a)-5(c) and 6(d)-6(f), for which only one supercritical cluster is present 
in the simulation box, the large and distinct jumps in R{t) are not clearly seen. 
Inspection of atomic distributions for these simulations shows that growth of this 
supercritical cluster is mainly realized via its fast diffusion among almost immobile 
individual copper atoms which are sometimes "swept" and absorbed by this mobile 
cluster. Some subcritical clusters containing several copper atoms are also observed 
in these simulations, and they seem to diffuse notably slower than the supercritical 
cluster, in a qualitative agreement with Fig. 9 in [lO|. 

The resulting acceleration diffusion parameter Ac for dilute Fe-Cu alloys is shown 
in Fig. 7. It is rather high: Ac > 200, having the same scale as the parameters of 
acceleration of precipitation in Fe-1.34%Cu alloys at T = 500°C due to the high 
mobility of copper clusters studied in Refs. [10|] and [21|]. Fig. 7 also shows that the 
concentration and temperature dependencies of this acceleration are rather sharp. 
Fig. 7(a) illustrates the notable increase of Ac with increase of copper concentration 
x; it can be explained by a probable enhancement of density of mobile copper 
clusters (both supercritical and subcritical) under increase of the copper content in 
an alloy. Fig. 7(b) shows a significant decrease of Ac under elevating temperature T. 
This also seems natural, as elevating T should lead to the weakening of the copper- 
vacancy binding proportional to the Mayer functions [exp (/3£'^q^) — 1] with E^q^ 
from (|4|) [18]. Hence the strong trapping of vacancies by copper clusters (which is 
the physical origin of their high mobility [l3]) should weaken. 



6. Comparison of nucleation barriers found using the STN calculations and 
the KMC simulations 

Estimates of the nucleation rate Jkmc described in Sec. 3 and illustrated by Figs. 
2 and 3, combined with the estimates of the effective diffusivity -Dg// discussed in 
Sec. 5 and illustrated by Figs. 5-7, enable us to find the reduced nucleation barriers 
f3Fc in Eq. ()18p basing on the kinetic Monte Carlo simulations. This equation (llSp 
also shows that for the large (3Fc considered which vary with x and T very sharply, 
possible errors of a relative order of unity in our estimates of parameters Jq, Jkmc 
and Dgff make no significant effect on the /3-Fc values obtained. 

In table 4 we present the (/3-Fc)kmc values estimated for all simulations shown 
in Figs. 2 and 3, together with the analogous (/3F^ )rtn va lues calculated in Sec. 2 
with the use of the statistical theory of nucleation jl2l . Il3l | . In comparison of these 
KMC and STN results we should remember that the STN assumes the size and 
the nucleation barrier of the embryo to be large: A'^c ^ 1) I^Fc ^ 1- Therefore, 
the accuracy of this theory and its agreement with the KMC simulations should 
improve when the size of the embryo increases. The results presented in tables 1, 2 
and 4 agree with these considerations. We see that at low concentrations x < 0.2% 
when both Nc and l3Fc are sufficiently large: Nc > 12, /SFc > 20, the agreement 
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between KMC and STN results is virtually perfect: for all nine simulations with 
X < 0.2% shown in table 4, differences between the (/3Fc)kmc and the (/3-Fc)stn 
values have the order of a percent. At the same time, at higher x equal to 0.3 or 
0.4%, the Nc and f3Fc values notably decrease, and differences between (/3Fc)kmc 
and {/3Fc)sTN increase, which may reflect the lowering of accuracy of the STN. 

Therefore, the results presented in table 4 enable us to make an important con- 
clusion that at high values /3-Fc ^ 20, the STN-based calculations of nucleation 
barriers are highly reliable. At the same time, as mentioned in Sec. 2, the homo- 
geneous precipitation limit in dilute Fe-Cu alloys corresponds just to the values 
l3Fc > 20. Therefore, the nucleation barriers near the homogeneous precipitation 
limits in these alloys can be reliably calculated using the STN. 



7. Temperature dependencies of precipitation limits in dilute iron-copper 
alloys 

Calculations of nucleation barriers described in Sees. 2 and 6 can be used for 
tentative estimates of the homogeneous precipitation limit x^°"^{T) at different 
temperatures T. As mentioned in Sec. 2, it seems natural to suggest that the 
position of this limit is mainly determined by the value of the activation factor 
ex.p{—(3Fc), that is, by the reduced nucleation barrier (3Fc which at x = xf°"^ takes 
a certain high value {(3Fc)i. It seems also natural to assume that these precipitation 
limits at different temperatures T correspond to the similar values of the activation 
factor exp(— /3Fc). It implies the following relation: 

(/3Fe)z = /3Fe [xf ""^ (T) , r] ^ Chom (37) 

where Chom is a constant independent of temperature, and thus it can be estimated 
using the experimental value xJ*°™'(290°C) = x^^^. For example, using estimates 
of x^^^ mentioned in Sec. 2, x^^^ ~ 0.08 — 0.09% from 3 or x^^^ ~ 0.15% from 0], 
we obtain: Chom — 28 or Chom — 23. Then equation (j37|) with the function Fc{x, T) 
calculated using the STN enable us to find the xf°™(T) value. A high accuracy 
of the STN calculations of nucleation barriers combined with the above-mentioned 
physical considerations allow us to expect that such "semi-empirical" estimates of 
the homogeneous precipitation limit for temperatures T of practical interest can 
be sufficiently reliable. 

One can also try to extend this approach to the case of the heterogeneous pre- 
cipitation of copper on clusters of irradiation-induced point defects, vacancies and 
self-interstitial atoms (that is, to the case of neutron or ion irradiation to 
estimate the appropriate heterogeneous precipitation limit xj^^*, which is more in- 
teresting for applications. Even though kinetic paths of the heterogeneous and the 
homogeneous precipitation differ from each other (see, e. g., [s^]), the main physical 
origin for suppressing precipitation at low solute concentrations x seems to be the 
same for both processes. In both cases, formation of a precipitate needs overcoming 
the energetic barrier due to the surface loss in the free energy F which is not com- 
pensated by the volume gain in F until the precipitate volume becomes sufficiently 
large. Differences between these two processes are usually related mainly to some 
geometrical factors which depend on the structure of the heterogeneity [32] • 

Therefore, one may expect that the heterogeneous precipitation limit x^^* can 
also be estimated from a phenomenological relation similar to ()37p . but the "lim- 
iting" value Fci in this relation should be higher as the nucleation barrier for the 
heterogeneous nucleation, Fj^'^^, is lower than that for the homogeneous nucleation, 
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phom_ However, if we assume that the difference between F^*^* and F^°^ has mainly 
a geometrical origin (which is the case for the simplest models of heterogeneous 
nucleation [13]), then the ratio F^""^ / F^""^ can be supposed to weakly vary with 
temperature. Then the temperature dependence of the heterogeneous precipitation 
hmit xf^*(r) can be estimated from a semi-empirical relation similar to ([31 



x^^\T),T ^Chet (38) 



where Fc{x, T) is again the STN calculated nucleation barrier for the homoge- 
neous nucleation, while the constant Ch^t is estimated using the experimental value 
xJ''^*(290°C) = x^g^p. For numerical estimates we use the value x^^^p ~ 0.06% 
given m Refs. and 0, which yields: Chet — 35. 

In figure 8 we show positions of the homogeneous and the heterogenous pre- 
cipitation limits in the (x,T) plane, TI^°"^{x) and T;'^^*(x), calculated for dilute 
Fe-Cu alloys using Eqs. (f37l) and (p8]l . For comparison, in Fig. 8 we also show 
the binodal curve Tb{x) (the copper solubility limit) estimated in from exper- 
imental data and ab initio calculations. The results in Fig. 8 are presented for the 
temperature interval between 290 and 450° C used for annealing of RP Vs [l, 0] . 
The value Ctom = 27.6 used in Fig. 8 corresponds to xJ"""(290°C) = 0.088% for 
which Radiguet et al. [6] did not observe any homogeneous precipitation, while 
the constant C^et = 35 used in Fig. 8 corresponds to the estimate xj**^* = 0.06% 
given by Miller et al. @, 0]- Our calculations also show that a slight variation of 
the constant Chom or Chet leads to an "almost rigid" shift of the curve Tj^°'^{x) 
or TI^^^{x) in Fig. 7. For example, when this variation of Chom leads to the shift 
,5a;f°'"=0.05% to the right at 290°C, the analogous shift at 450°C is almost the 
same: <5xJ"""(450°C) ~ 0.06%. 

The most interesting general feature seen in Fig. 8 seems to be a rather sharp tem- 
perature dependence of the precipitation limits, particularly for the heterogeneous 
precipitation. For example, the x^^^ value at T = 450° C exceeds that at T = 290° C 
by almost five times, while the x^°^ value increases for this temperature interval 
by about 3.5 times. These qualitative conclusions can be useful, in particular, for 
the interpretation of microstructural observations related to the annealing of RPVs 
[3]. Note also that the variations of both precipitation limits, x^^^{T) and x^°^{T), 
with temperature for the temperature interval considered are much stronger than 
those for the solubility limit Xb{T). 

To conclude this section, we again note that our equations (f3711 and ([38]) for 
temperature dependencies x|^^*(r) and x^°"^{T) are of course tentative and have 
not been proved formally. At the same time, physical considerations about the 
dominant role of the activation factor exp {—f3Fc) in suppressing precipitation at 
low X used in our derivation seem to be rather plausible, particularly for the ho- 
mogeneous precipitation. In the derivation of Eq. (I38p for the heterogeneous pre- 
cipitation, we also used an assumption about a weak temperature dependence of 
the ratio F^^^/F^°^ for the temperature interval considered which may seem to 
be less evident. Therefore, the accuracy of Eq. dMl) for x^''\T) can, generally, be 
lower than the accuracy of Eq. dniD for xf°'^(r). However, one may expect that 
the main qualitative features of both dependencies, Tj^°^{x) and TI^'^^{x), shown 
in Fig. 8 are correctly described by the simple model used. 
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8. Conclusions 

Let us summarize the main results of this work. We study kinetics of homogeneous 
nucleation and growth of copper precipitates under electron irradiation of iron- 
copper alloys at low concentrations x = 0.06 — 0.4 at.% and temperatures T = 
290 — 450°C used for service of a number of nuclear reactor pre ssure vessels [H-0]. 



The earlier-described kinetic Monte Carlo (KMC) modeling 10(] and the statistical 
theory of nucleation (STN) (l^ . [l^ are used. The ab initio model of interatomic 
interactions which describes well the available data about precipitation in Fe-Cu 
alloys at different x and T [lol . 11] is used for both the KMC simulations and 
the STN calculations. The conventional assumption [2] about the similarity of 
mechanisms of precipitation under electron irradiation and under thermal aging is 
also adopted. Then precipitation under electron irradiation can be described by the 
KMC codes developed by Soisson and Fu [13] for studies of precipitation during 
thermal aging but with the acceleration of kinetics under irradiation characterized 
by the acceleration factor Airr in Eqs. ([2]) and ([3]). We estimate this acceleration 
factor for dilute iron-copper alloys considered, and our estimate reasonably agrees 
with the available experimental data Q]- 

Our STN-based calculations of the nucleation barrier Fc in the Zeldovich-Volmer 
relation ([1]) for the nucleation rate J show that this nucleation barrier varies with 
X and T very sharply. Thus the concentration and temperature dependencies of 
the nucleation rate are mainly determined by the variations of the activation fac- 
tor exp(— /?Fc) with x or T. We also found that at T = 290°C, the interval of 
concentrations x for which this activation factor starts to fall off very rapidly just 
corresponds to the interval of positions of the homogeneous precipitation limit xj""" 
estimated in experiments [1, 0, 0] • 

Our KMC simulations for the dilute alloys considered describe mainly the very 
initial stage of precipitation, that of the steady-state nucleation characterized by 
the constant nucleation rate J. We also use these KMC simulations to study kinetics 
of growth of a supercritical embryo and confirm the validity of the phenomenolog- 
ical equation of the classical theory of nucleation which describes this growth [141 ]. 
Our simulations also enable us to estimate the effective diffusivity -De// which en- 
ters into the STN expression for the nucleation rate J. The -De// values are found 
to exceed the dilute alloy values Dq{T) by two-three orders of magnitude. This 
effective acceleration of diffusion seems to be mainly due to the high mobility of 
small copper clusters found earlier for the precipitation in Fe-1.34%Cu alloy at 
T = 500°C [13, Eli. 



The KMC estimates of the nucleation rate J and the effective diffusivity -De// 
described above enable us to determine the nucleation barrier Fc using Eq. (jlSp 
based on the KMC simulations. The resulting values of the reduced nucleation 
barrier (/3Fc)kmc ai'e compared with the (/3-Fc)stn values calculated using the 



statistical theory of nucleation [Ij, ll3|]. We find that for the dilute alloys with 
X < 0.2%, the STN and the KMC results for the nucleation barrier Fc coincide 
within about a percent. It seems to confirm a high reliability of the STN for this 
problem. 

Making a plausible assumption that the position of the homogeneous nucle- 
ation limit x^""^ is mainly determined by the value of the activation factor 
exp [— /3Fc(a;, T)] which at x = x^°'^ takes an approximately same value exp {—f3Fc)i 
for all temperatures T considered, we use our STN-based calculations of reduced 
nucleation barriers jSFc to estimate the temperature dependence x^°"^{T). For tem- 
peratures between 290 and 450°C, this dependence is found to be rather sharp, 
much sharper than that for the solubility limit XgoiiT). Making an additional as- 
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sumption about a weak temperature dependence of the ratio of nucleation barriers 
for the heterogeneous and the homogeneous nucleation, Fj^^^/F^'^, we also esti- 
mate the temperature dependence of the heterogeneous precipitation limit x^'^*{T) 
which corresponds to the neutron or ion irradiations. The dependence x[^^*(T) is 
found to be still more sharp than that for the homogeneous precipitation. In spite of 
the evidently tentative character of these estimates, one can expect that the main 
qualitative features of both dependencies, x^°^{T) and xf'^*(T), are described by 
this model properly. 
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Table 1. Some characteristics of critical embryos in alloys Fei-^Cui; at T = 290° C 
calculated using the "pair-cluster-with-fluctuations" approximation of STN [l^] 
and Cu-Cu interactions given by Eqs. A'^c is the number of copper atoms within 
the embryo, and Rc is its effective radius defined by Eq. ([8]); Af^o and Ai^i are the 
nucleation barrier values found in the zeroth and the first order in fluctuations, as 
explained in the text, and = AJIq + Aili is the full nucleation barrier. 
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Table 3. Values of various factors in the expression ([9|) for the prefactor Jq in Eq. ([T]) 
for some Fei-^Cua; alloys; 'jnn is the effective exchange rate between neighboring 
copper and iron atoms given by Eq. (jl4p . and Va = ao/2 is volume per atom in the 
BCC lattice of iron. 
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Table 4. Values of the reduced nucleation barrier in Fe-Cu alloys estimated in our 
KMC simulations, (/3Fc)kmc! and those calculated using the STN, (/3-Fc)stNi at 
different temperatures T, concentrations x, and simulation sizes L„ in (jl5p . 
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Figure 1. (color online) Concentration profile Ac(ri) = c(ri) — c in the critical embryo; c{ri)=Ci is the 
mean occupation of site i by a copper atom (local concentration) at the distance from the center of 
embryo, and c = x is the copper concentration far from the embryo. Solid or dashed lines are drawn 
to guide the eye. Curves A, B, C and D (red, green, blue and grey online) correspond to the following 
temperatures T (in °C) and concentrations x (in at.%). A: T = 290, x = 0.06; B: T = 290, x = 0.2; C: 
T = 290, X = 0.3; and D: T = 390, x = 0.2. 

Figure 2. (color online) Temporal evolution of the precipitate density dp = Np/Vs observed in our KMC 
simulations with the final number of precipitates Np > 3. Simulation time on abscissa axis is given in the 
"thermal aging" values t^h related to the observed time tirr by Eq. l(3]|. Each of frames (a) -(f) shows 
simulation for the following concentration x (in %), temperature T (in °C) and the simulation size L in Ln 
given by Eq. ^Ej/: (a) 2:=0.172, r=290, L = L2; (b) a:=0.172, T=290, L = L3; (c) x=0.2, T=290, L = Li; 
(d) x=0.2, T=290, L = L3; (e) a;=0.3, T=290, L = Lr, and (f) x=0.3, T=290, L = L3. For simulations 
shown in frames (e) and (f) we also present the mean precipitate radius Rm obtained by averaging of the 
left-hand side of Eq. I I27I1 over all precipitates within simulation box. Dashed lines (blue online) show our 
estimates of the coefficient J = Jkmc in Eq. I|16|l . 

Figure 3. (color online) The same as in Fig. 2 but for the following simulations: (a) x=OA, T=290, L = L\ \ 
(b) x=0.2, r=305, L = Lr, (c) a:=0.2, T=305, L = L2; (d) a;=0.2, r=305, L = L3; (e) x=0.2, T=320, 
L = Li; and (f) x=0.2, T=320, L = L3. 

Figure 4. Evolution of point defect atomic fractions in pure iron, during an irradiation at 2 X 10~^ dpa.s"^ 
and T = 290°C, with a dislocation densities (a) p^j = 10* cm^^ and (b) = 10^^ cm~^. Solid line (red 
online): vacancies; dashed line (green online): interstitial atoms. 

Figure 5. (color online) Solid lines (red online) show the dependencies R = RKMci^th) with R defined 
by Eq. l|27|l observed in our KMC simulations for temperature T = 290°C. For simulations with A'p = 1 
(those for x < 0.172%), this dependence corresponds to the single precipitate within simulation box; and 
for simulations with Np > 1 (x > 0.172%), to the biggest precipitate within simulation box. Dashed line 
(blue online) shows the dependence R{^) given by Eq. I I31II . Each of frames (a) -(f) corresponds to the 
following values of the concentration x (in %), the simulation size L (in Ln), and the acceleration diffusion 
parameter A^ in Eq. (a) x=0.088, L = Li, Ac = 334; (b) 3;=0.116, L = Li, Ac = 396; (c) 2:=0.144, 

L = Li, Ac = 450; (d) x=0.172, L = L3, Ac = 570; (e) a;=0.2, L = Li, Ac = 463; and (f) x=0.3, L = L3, 
Ac = 666. 

Figure 6. (color online) The same as in Fig. 5 but for the concentration x = 0.2%. Simulations for 
T < 320°C correspond to Np > 1, and those for T > 320°C, to Np = 1. Each of frames (a) -(f) 
corresponds to the following values of temperature T (in °C), the simulation size L (in L„), and the 
acceleration difl'usion parameter Ac. (a) T=305, L = Li, Ac = 430; (b) T=305, L = Lg, A = 350; (c) 
r=320, L = Li, Ac = 330; (d) r=335, L = Li, Ac = 200; (e) r=380, L = L-^, Ac = 230; and (f) T=390, 
L = Li, Ac = 200. 

Figure 7. (color online) (a) Values of the acceleration diffusion parameter Ac in Eq. l|34|l observed in our 
KMC simulations for Fei-^^Cu^; alloys at temperature T = 290°C and different concentrations x. Circles, 
squares and triangles (green, blue and red online) correspond to the simulation size L equal to Li , L2 and 
L3, respectively, (b) The same as in (a) but at x=0.2% and different temperatures T. 

Figure 8. (color online) Solid line (red online) shows the homogeneous precipitation limit for dilute Fe-Cu 
alloys, Tl^°"^{x) or x^°"^(T), estimated using Eq. JSTjl with Chom = 27.6. Dashed line (green online) shows 
the heterogeneous precipitation limit for these alloys, T^^^{x) or x[^'^^(T), estimated using Eq. (jsHl with 
Chet = 35. Chain line (blue online) shows the binodal Ti,(x) as estimated in [To| . 
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